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A MODEL FOR TAX EVASION WITH SOME REALISTIC PROPERTIES 


RICHARD VALE 


Abstract. We present a discrete-time dynamic model of income tax evasion. The model is solved exactly in the 
case of a single taxpayer and shown to have some realistic properties, including avoiding the Yitzhaki paradox. The 
extension to an agent-based model with a network of taxpayers is also investigated. 


1. Introduction 

1.1. There are various reasons for building a mathematical model of tax evasion. It is interesting to create 
models which exhibit some features of the real world, as this might help us to understand the mechanisms 
underlying taxpayer behaviour. It is also useful to have models available for testing the possible effects of new 
compliance strategies such as changes in tax rates. 

1.2. We present a model of income tax which combines new ideas with ideas from existing models. The model is 
dynamic (evolves in time) and quite simple, having four parameters in the one-taxpayer case. These parameters 
are the tax rate, the probability of audit, the penalty rate and the savings rate (see Section]^ for the details.) The 
new model is interesting because it has two realistic features. 


1.3. Firstly, many existing models predict that there should be widespread non-compliance, but in fact many 
people pay tax. Some authors explain this by introducing a notion of tax morality or proposing that there is a 
“social reward” for paying tax which should be included in calculations. Since taxes pay for infrastructure, it 
makes sense that taxpayers would be motivated by social considerations. Nevertheless, the present paper aims to 
show that widespread compliance can arise in a population of purely selfish taxpayers. 

1.4. Secondly, another strange aspect of tax evasion models is that many of them predict that increasing the tax 
rate should increase compliance. This does not seem to agree with real life; when taxes are very high, they are 
viewed as unfair, and people are more reluctant to pay them. In the extreme case, everybody would be expected to 
evade a tax over 100/l[^as this is equivalent to confiscation of all income. The problem of models predicting that 
higher taxes should result in higher compliance is sometimes called the Yitzhaki Problem or Yitzhaki Paradox. 

1.5. The model presented in this paper predicts that higher taxes should result in lower compliance, provided 
that the tax rate exceeds a certain critical value, and a 100% tax will result in total non-compliance. Again, the 
emotions of the taxpayers do not need to be taken into account. The relationship between tax rate and 
non-compliance is derived from economic considerations and taxpayers seek only to maximise their prohts. 


1.6. In the simplest terms, the reason why the new model predicts widespread compliance is because the tax 
authority is assumed to collect all unpaid taxes rather than just taxes from the most recent period. Some authors 
have previously considered models like this, although Seibold and Pickhardt (who call this feature back-auditing 
0) comment that it has been ‘largely neglected in the literature’. The reason why the new model predicts higher 
non-compliance at higher tax rates is that people will not have enough money to live on if the tax rate is too high. 


1.7. References. The most famous model of tax evasion is the 1972 Allingham and Sandmo model [T]. The 
Yitzhaki Paradox originates in Yitzhaki’s 1974 paper [6]. The question of the effect of tax rates on tax compliance 
in the real world has been much studied, and a survey of models, empirical studies and experimental results is [2]. 
The model introduced in the present paper is related to a model of Klepper, Nagin and Spurr [1] in which 
taxpayers make a decision whether to evade based on their savings rate. However, the model of [3] only has two 
periods {t = 0,1.) The original aim of the current paper was to develop a model for a single taxpayer which can 
be used as an agent in an agent-based model of tax evasion on a network. This was inspired by the agent-based 
model of Korobow, Axtell and Johnson [3]. The model using a toroidal grid in Section 5.7 is closely based on the 
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^We consider only simple models in which a taxpayer’s entire income is taxed at a fixed rate. It is not clear that taxpayers with a 
marginal tax rate of 100% would be expected to evade. 
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Ising model of statistical mechanics. The idea of using the Ising model to study tax evasion is due to Zaklan, 
Westerhoff and Stauffer [7]. 


2. A Model of a Single Taxpayer 

2.1. In this section we present a model for a single taxpayer. The model turns out to be simple enough to have 
an explicit solution. First we describe the model in words. There are four parameters: the tax rate r, the savings 
rate k, the probability of audit p and the penalty rate A. We assume 0 < T,k,p < 1 and A > 1. 

2.2. Verbal Description. The taxpayer has an income of $1 and can choose whether to pay the tax of $r. If the 
taxpayer pays the tax, the taxpayer is said to comply and receives $(I — t). Otherwise, the taxpayer evades and 
receives $1. The taxpayer saves a proportion k of this income and spends the rest. The taxpayer may then be 
audited. If audited, the taxpayer must pay back all the tax evaded so far, together with a penalty which is 
proportional to the amount of tax evaded. Because the taxpayer spends some of his or her income, it may not be 
possible to recover all of the evaded tax. In this case, the taxpayer pays as much as possible, and then continues 
with $0 at the next time step. 

2.3. The decision whether to evade is based on the profit made from evasion so far. If the profit made from 
evasion between time 0 and time t is positive then the taxpayer evades at time t, otherwise the taxpayer complies 
at time t. 


2.4. The model is a very simplified version of how a real-world tax authority operates. The tax authority, upon 
uncovering evidence of evasion, attempts to recover all the evaded tax, not just the tax from the most recent 
period. But this is not always possible, because the taxpayer might have spent some of the evaded tax. Real-world 
tax agencies might seize the taxpayer’s assets, but might still fall short if the evaded tax had been spent on 
perishable items such as meals, holidays or gambling. Interest is ignored in the model. In real life, the taxpayer 
might earn interest on savings, but this is balanced by the interest that must be paid on evaded tax. Of course, 
this is one of many assumptions that simplify the model. For example, we ignore utility and assume that tax is 
paid on the whole income at the same fixed rate. 


2.5. Detailed Description. We now give a more detailed mathematical description of the model. Time moves in 
discrete steps t = 0,1, 2,.... At time t, the taxpayer may or may not evade (the condition for evasion will be given 
below.) If the taxpayer evades, the income received is 1. If the taxpayer does not evade, the income received is 
1 — r. A proportion k of the income is retained. Let f{t) denote the taxpayer’s fortune at time t. Then 


/(^ + 1 ) 


fit) + k 

fit) + kil - t) 


if taxpayer evades at time t 
otherwise. 


Denote by pf(t) the profit made from evasion up to and including time t. Then 


pf(i + 1) 


pf(t) -br 
pf(t) 


if taxpayer evades at time t 
otherwise. 


Denote by n(t) the number of times the taxpayer has evaded since the last audit. Then 


n{t -b 1) = 


n(t) -b 1 
n{t) 


if taxpayer evades at time t 
otherwise. 


Suppose the taxpayer is audited at time t with probability p. Audits are independent of one another. If the 
taxpayer is audited then the taxpayer must pay ATn(t -b 1). If fit -b 1) < ATn(t -b 1) then the taxpayer must pay 
the entire fortune fit + 1) instead. Thus, if audited, /, pf and n are updated according to the rules 


fit + 1) ^ fit + 1) - min(/(t-b 1), ATn(t-b 1)) 
pf(t + 1) ^ pf(t + 1) - niin(/(t -b 1), At n(t -b 1)) 
n(t -b 1) •«— 0 


It remains to give the condition for evasion. We assume that the taxpayer chooses to evade if pf(t) > 0. In other 
words, if the profit made from evasion so far is positive, the taxpayer chooses to evade because evasion looks like a 
good bet. We also assume /(O) = n(0) = 0 and pf(0) > 0 (otherwise the taxpayer will never evade.) 
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2.6. Difference Equations. The model can be described even more formally via the following system of 
equations 

At ~ Bernoulli (p) 

fit + 1) = fit) + fc(l - "t) + fcT(5pf(t)>o - min (/(t) + fc(l - r) + hrSp^t) > 0 , AT(n(t) + (5pf(t)>o)) At 
pf(t + 1 ) = pf(t) + T - min (/(t) + k{l - r) + kT6pi(^t)>o, AT(n(t) + (5pf(t)>o)) At 
n(t + 1 ) = (n(t) + (5pf(t)>o)(l — At) 

where 5pf(t)>o is 1 if pf(t) > 0 and 0 if pf(t) < 0, and /(O) = n(0) = 0 and pf(0) > 0. 


3. Solution of the One-Taxpayer Model 

3.1. In this section we describe the dynamics of the one-taxpayer model. If pf(t) < 0 at some time t, then the 
taxpayer never evades again, so we restrict to the case pf(t) > 0. In this case, pf(t) cannot decrease until the 
taxpayer is audited, so between audits the taxpayer always evades. 


3.2. Suppose the taxpayer is audited at time Tq and at time Ti > Tq and there are no audits in between. We 
consider the total profit made from evasion 

pf(ri) - pf(To) 


between Tq and Ti. 


3.3. In the first case, suppose At < k. Then f{t) > ATn(t) for all t. This follows by induction on t. It is true for 
t = 0. If f{t) > ATn(t) and At = 0 then + 1) = fit) + k > Ar n(t) -|- Ar = ATn(t -1-1). If fit) > Ar n(t) and 
At = 1 then n(t -|- 1) = 0 and so f{t -I- 1) > XTn{t + 1) in this case too. Therefore, in case At < k, we have 

fit) > ATn(t) for all t and 

pf(Ti) - pf(To) = t(Ti - To) - min(f(Ti) + k, AT(n(Ti) -P 1)) 

= T(ri-To)-AT(n(Ti) + l) 

but n(ri) -I- 1 = Ti — To because the taxpayer evades at times Tq -I- 1, Tq -f 2,..., Ti and so 

pf(Ti) - pf(To) = t(1 - A)(Ti - To) 

3.4. In the second case, suppose At > k. Then f{t) < ATii(t) for all t. This follows by induction on t. It is true 
for t = 0. If f{t) < ATn(t) and At = 0 then fit+ 1) = /(t) -I-A: < ATn(t) -|- At = ATn(t-|- 1). If /(t) < ATn(t) and 
At = 1 then f{t -f 1) = 0 and so 0 = f{t -I- 1) < ATn{t + 1) in this case too. Therefore, in case At > k we have 
fit) < At nit) for all t and 

pf(Ti) - pf(To) = t(Ti - To) - min(/(Ti) -f k, AT(n(Ti) -P 1)) 

= T(Ti-To)-(/(Ti) + fc) 

but /(Ti) + k = fc(Ti — To) because /(To) = 0 in this case (the taxpayer is audited at time To by assumption and 
loses their entire fortune because /(To — 1) -I- fc < AT(n(To — 1) + 1)) and the taxpayer evades at times 
To -|- 1, To -|- 2,..., Ti, and so 

pf(Ti)-pf(To) = (T-A:)(Ti-To) 


3.5. Although To and Ti are random variables, we see that in all cases the quantity 


pf(Ti) - pf(To) 
Ti - To 


= T — min(fc, At) = 



T > k/A 
T < k/A 


is constant. This quantity represents the average change in pf between audits. Equivalently, it represents the 
average profit from evasion. We will use it extensively in this paper and denote it using the derivative-like symbol 

Apf 
AT ■ 
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tau = 0.02 tau = 0.39 tau = 0.42 



Figure 1. Plots of pf(t) versus t from simulations from the one-taxpayer model with k = 0.4, 
p = 0.01, pf(0) = 5, A = 1.5 and r = 0.02,0.39,0.42. The dotted lines are y{t) = pf(0) -|-t(1 — X)t 
(left) and y(t) = pf(0) -f (u — k)t (middle and right). 


3.6. We see that if r > fc then pf is an increasing function from one audit to the next, and therefore the taxpayer 
always evades. If r < fc then pf will eventually become negative and the taxpayer becomes compliant. In this case, 
we can calculate the expected amount of time which elapses until the taxpayer becomes compliant. Let T^omp be 
the time at which pf first becomes negative. Each audit reduces pf by x AT where AT is the time elapsed 
since the last audit. So pf will become negative as soon as there is an audit after time — pf(0)/^^ has elapsed. 
The expected time until the next audit is a geometric random variable with a mean of 1/p and therefore 

. pf(o) 

E[Tcomp] 


3 ] = -pf(0) 


Apf\ 


AT ) 


1 


-b - = 


- J 


pf(0) 

. '^(>‘- 1 ) 


T > fc/A 
r < fc/A 


(3.1) 


In particular, the expected time until the taxpayer becomes compliant is a linear function of — 

3.7. Examples. The calculations of Section]^ can easily be tested by implementing the equations of Section 2.6 
on a computer. Plots of pf(t) against t from three example runs are shown in Figure]^ Here fc = 0.4 and A = 1.5. 
The first two graphs show results for t = 0.02 and r = 0.39. These taxpayers eventually become compliant but at 
different rates, as shown by the dotted lines with slope The third taxpayer has r = 0.42. This tax rate is too 
high for the taxpayer ever to become compliant. 

3.8. Measuring Evasion. When t < fc, the taxpayer always becomes compliant eventually. However, Equation 
3.1 shows that the expected time taken to reach compliance (which is the same as the number of times the 
taxpayer is expected to evade until becoming compliant) is a linear function of —(Apf /AT)~^ with positive 
coefficients. It is therefore reasonable to take —(Apf /AT)~^ as a measure of non-compliance. It is likewise 
reasonable to take as a measure of non-compliance because it measures how profitable evasion is for the 
taxpayer on average. 


3.9. The one-taxpayer model is interesting because it exhibits phenomena which can be seen in the real world 
but not in many models of tax evasion. A graph of against t for fixed values of fc and A is shown in Figure 
Also sketched in Figure[^is a graph of —(Apf /AT)~^ against r. It will be useful to compare the U-shape of this 
graph with experimental results from Figures and 

3.10. Discussion. If r > fc/A, increasing the tax rate causes compliance to decrease. This fact does not require 
any arguments about the morality of taxpayers. It is a consequence of the cost to the taxpayer of non-compliance 
being bounded below, which means that evasion becomes more profitable to the taxpayer when the tax rate is 
higher. 

3.11. When T ^ fc, the taxpayer becomes compliant. When t ^ fc/A, increasing the tax rate causes 
non-compliance to decrease. Non-compliance is minimised when r = fc/A. This optimal tax rate increases with the 
rate of saving and decreases when the penalties for evasion become harsher. 
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minimised when r = fc/A. Right panel: graph of —(Apf/AT) The shape of this graph should 
be compared with the experimental results in Figure]^ and Figure]^ 


3.12. Progressive taxation. The result that a tax rate of r = fc/A minimises non-compliance might help to 
justify progressive taxation (which is a form of taxation in which income over a certain threshold is taxed at a 
higher rate.) All other things being equal, it seems reasonable that a higher proportion of income will be saved if 
it is above what is needed for survival. This corresponds to a higher fc and therefore a higher fc/A. This argument 
is similar to a common argument in favour of progressive taxation based on marginal utility but has a different 
motivation; the model suggests that taxes should be progressive if the tax authority wishes to minimise 
non-compliance. 


4. A Network of Taxpayers 

4.1. In the real world, there is more than one taxpayer. One of the original motives for creating the model of 
Section was to use it as a model for a single taxpayer operating within a network of taxpayers who share 
information. The idea of investigating such a model was inspired by the work of Zaklan, Westerhoff and Stauffer 
[7] and the agent-based model of Korobow, Axtell and Johnson [3]. In this section, we construct a model for a 
network of taxpayers and present some experimental results. 

4.2. We assume that N taxpayers are connected in an undirected network. Let A^y = I if a: and y are connected 
and Axy = 0 otherwise. Assume that A^x = 0 for all x (there are no self-loops.) Suppose the tax rate t , savings 
rate fc and probability of audit p are fixed quantities in (0,1), and there is a fixed penalty rate A > I. 

4.3. For a taxpayer x, define the neighbourhood 

N{x) = {a;} U {y : Axy = 1} 

to be the set of taxpayers consisting of x and the neighbours of x. 

4.4. Time proceeds in discrete steps t = 0,1, 2,.... Let f{x, t) be the fortune of taxpayer x at time t, pf(a:, t) be 
the profit made by x from evasion up to and including time t, and n(a:, t) be the number of times x has evaded 
since the last audit. Assume that taxpayers are audited independently of one another and that the probability 
that a given taxpayer x is audited at a given time t is p. Assume /(x, 0) = n(x, 0) = 0 for all x and that /, pf and 
n evolve for each taxpayer exactly as in the one-taxpayer model of Section except for the criterion for choosing 
to evade. 
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node 

1 (centre) 

2 

3 

4 

5 

6 

7 

8 

9 

10 

mean 

420 

410 

394 

434 

435 

408 

415 

412 

408 

425 

sd 

153 

239 

240 

248 

242 

237 

236 

254 

231 

237 


Table 5.1. Time until compliance for ten taxpayers in a star-shaped network for 50 simulations 
of the model. See Section [5^ for details. 


4.5. In the new model, the network structure enters through the rule which a taxpayer uses to decide whether to 
evade. A taxpayer is assumed to be aware of their own history and the histories of all of their neighbours and the 
decision to evade is based on their best guess at the expected proht from evasion given this information. Let 
ntotai( 2 ;, t) be the total number of times that x evades between time 0 and time t. A taxpayer x chooses to evade if 

ntotai(y,t) > 0. (4.1) 

y£N{x) y£N{x) 

which is the same as 

pf(?/,t)>0. (4.2) 

y^N{x) 


5. Analysis of the network-of-taxpayers model 

5.1. The model for a network of taxpayers is no more difficult to implement than the one-taxpayer model. 
Simulations show that a network of taxpayers behaves in a similar way to a single taxpayer from a qualitative 
point of view. However, it seems to be much more difficult to calculate exactly what will happen for a network in 
the manner of Section The questions of interest are: Do the taxpayers eventually become compliant? How long 
does it take? Does the time until compliance vary from one taxpayer to another? 

5.2. For a taxpayer x, pf(a:,t) evolves in much the same way as in Section]^ In particular, we need only consider 
the case < 0. But it is no longer guaranteed that pf(a:,t) will become negative as t ^ oo; the taxpayer may 
become compliant before this happens. Nevertheless, if taxpayer x is non-compliant, then pf(ai,t) will decrease 
when X is audited. Therefore, the taxpayers must eventually become compliant. 


5.3. We make only a crude estimate of the time until compliance. Let X be the set of taxpayers and A (Z X. 
(The case A = N{x) is of particular interest.) Define 

pf(Ai) = Y 

x^A 


Ignoring the randomness in the audits and assuming that every taxpayer evades at every time step, we have 


pf(x,t) « pf(a:,0) -I- 



t 


and so 

pf(A, t)psY Pf(^> 0) + 1^1 (^ 

x^A ^ ^ 


which becomes negative at time 




ji[ExGAPf(a^.O) 



(5.1) 


Therefore, if all the pf(a:, 0) are assumed to be equal, every subset of the taxpayers will become compliant at the 
same time. In particular, all the taxpayers are expected to become compliant at the same time. The time taken to 
become compliant does not depend on the network topology. 


5.4. Example. The model was simulated for a star-shaped network with ten nodes labelled with the integers 1 to 
10. Node 1 was connected to all other nodes and there were no other connections. Taking pf(x, 0) = 1 for all x 
and k = 0.4, r = 0.3, A = 1.5, p = 0.01 the results of 50 simulations are shown in Table [O] For each node x, the 
time in the table is the number of the last iteration on which x evaded. 
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Figure 3. Average number of evaders versus tax rate r for a star-shaped network with 10 taxpayers 
run for 1000 iterations. The line is at r = fc/A, the tax rate which minimised non-compliance in 
the one-taxpayer model. Compare Figure 


5.5. There is no evidence of a difference between nodes in the times until the nodes become compliant. However, 
Equation 5.1 gives an estimate of l/(fc — r) = 10 for the time for a node to become compliant, which is clearly 
wrong. This is because Equation |5. 1 1 ignores the audit probability p, which influences the time until compliance. 
Giving an accurate estimate of the time until compliance seems to be a more difficult problem than in Section 
especially when the pf(x, 0) are allowed to take arbitrary values. 


5.6. Figure shows the average number of evaders plotted against the tax rate for the star-shaped network with 
10 vertices and k = 0.4,p = 0.01, A = 1.5. The vertical line is at the value t = k/X. Comparison with Figure]^ 
shows that there is a wide range of parameters in the network case for which non-compliance is approximately 
minimised. All values of r between about 0.1 and 0.3 seem to be equally good. 


5.7. Example. Zaklan, Westerhoff and Stauffer [7] considered an Ising-like model in which taxpayers are 
connected in a square lattice. Korobow, Johnson and Axtell [3] considered a different model with a similar 
network topology but with diagonal connections. We simulated the model for 100 taxpayers arranged in a square 
lattice. The lattice is wrapped around to form a torus, so for example the taxpayer at (1,1) is adjacent to (1, 2), 
(2,1), (1,10) and (10,1). This ensures that every taxpayer is adjacent to 4 other taxpayers. As in the previous 
example, it is predicted that a tax rate of t = fc/A will minimise non-compliance. 

5.8. A plot of the average number of evaders against the tax rate r is shown in Figure]^ It can be seen that the 
shape is very similar to Figure which is further evidence to suggest that the topology of the network has no 
influence on compliance behaviour in this model. 


5.9. Varying the probability of evasion. We can make the taxpayer model even more like the Ising model by 
specifying a probability of evasion. The probability of taxpayer x evading at time t could be chosen to be an 
increasing function of J2y^N(x) By analogy with the Ising model, this could for example be taken to be 
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Figure 4. Average number of evaders versus tax rate r for a toroidal grid with 100 taxpayers run 
for 1000 iterations. The line is at r = fc/A, the tax rate which minimised non-compliance in the 
one-taxpayer model. Compare Figure and Figure 


where /3 > 0 is a parameter which plays the role of the reciprocal of temperature. As /3 —>■ oo, this reduces to 
(4.21. For smaller values of /3, there is more randomness in the decisions of the taxpayers and as t —>■ oo, the model 
approaches Bernoulli trials with a probability exp(/?^^) of evasion. Otherwise, this choice of evasion probability 
does not affect the qualitative behaviour of the model when the parameters are held constant, so we choose not to 
discuss it further. 


5.10. EfTect of heterogeneous k. The model is still unrealistic in many ways. For example, real taxpayers do 
not inhabit a torus, and the savings rate k undoubtedly varies over time and from person to person. There are 
many directions for further investigation. In this section we investigate what happens when each taxpayer is 
allowed to have an individual value for k. This follows the ideas of Korobow, Axtell and Johnson, who allowed 
various parameters to vary across taxpayers in their model. 

5.11. Under the conditions of Section [5?7l suppose that the savings rates of the 100 taxpayers are drawn 
independently from a Beta(2,3) distribution. This ensures that the mean savings rate is still about 0.4. The tax 
rate is chosen to be r = 0.3 and the probability of audit is p = 0.1. The model is run for 10000 iterations. 
Everything else is kept the same as in Section [STj We investigate the effect of varying k on compliance. We 
expect that individuals with higher k should be more compliant on average, but this might be tempered by 
network effects. 

5.12. It is of interest to see how many times a given taxpayer evades in the 10000 iterations. Figure shows a 
plot of a typical model run. The 100 taxpayers are plotted in a 10 x 10 grid and the darker squares indicate 
taxpayers who evade more often. The network is a torus, so the bottom row and top row are adjacent and the left 
column is adjacent to the right column. Note that the black and white regions in Figure]^ are almost contiguous, 
which is surprising because the savings rates of individual taxpayers were chosen at random. However, 
experiments with larger grids, for example 100 x 100, suggest that evaders tend to occur in small cross-shaped 
regions which are scattered at random throughout the network. 

5.13. The number of evasions for the 100 taxpayers is plotted against the value of k in the left panel of Figure 
There seems to be no obvious pattern. In the right panel of Figure the number of evasions is plotted against the 
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Figure 5. Total number of evasions for each taxpayer in a 10 x 10 toroidal grid run for 10000 
iterations. The value of k varies from taxpayer to taxpayer. The other parameters have the values 
A = 1.5,r = 0.3,p = 0.1. 
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Figure 6. Total number of evasions plotted against k (left panel) and average fc in a five-taxpayer 
neighbourhood (right panel) for each taxpayer in a 10 x 10 toroidal grid. The value of k varies from 
taxpayer to taxpayer. The other parameters have the values A = 1.5, r = 0.3, p = 0.1. The vertical 
lines are at k = t. 
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average value of k, fcavg = J2y£N{x) taxpayers x. It seems that taxpayers with fcavg(a:) < t evade 

every time, but there is no obvious pattern for taxpayers with k^^vgix) > t. Also, the majority of taxpayers either 
become highly compliant or highly evasive, with few in between the two extremes. 

6. Acknowledgements 
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7. Conclusion 

7.1. We have given a new model for tax compliance. The model for a single taxpayer is simple enough to be 
exactly solvable and has two realistic properties: it predicts that sensible tax rates will eventually lead to 
compliance and that there is a realistic relationship between changes in the tax rate and changes in the level of 
non-compliance. We propose that variants of this single-taxpayer model might be useful as agents in an 
agent-based tax evasion model. 
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